Rank-based model for weighted network with hierarchical organization and 

disassortative mixing 
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Motivated by a recently introduced network growth mechanism that rely on the ranking of node 
prestige measures [S. Fortunato et al, Phys. Rev. Lett. 96, 218701 (2006)], a rank-based model for 
weighted network evolution is studied. The evolution rule of the network is based on the ranking 
of node strength, which couples the topological growth and the weight dynamics. Both analytical 
solutions and numerical simulations show that the generated networks possess scale- free distributions 
of degree, strength, and weight in the whole region of the growth dynamics parameter (a > 0). We 
also characterize the clustering and correlation properties of this class of networks. It is showed that 
at Q = 1 a structural phase transition occurs, and for q > 1 the generated network simultaneously 
exhibits hierarchical organization and disassortative degree correlation, which is consistent with a 
wide range of biological networks. 

PACS numbers; 89.75.-k, 89.75.Hc 



A major source of the recent surge of interest in com- 
plex networks has been the discovery that a large num- 
ber of real-world networks have a power-law degree dis- 
tributions, so called scale-free networks [1-4]. Due to 
the peculiar structural features and the critical dynami- 
cal processes taking place on them [5-8], there has been 
a tremendous number of works modeling networks with 
scale-free properties. The previous models of complex 
networks always incorporate the preferential attachment 
[4], which may results in scale- free properties. That is, a 
newly added node is connected to preexisting one with a 
probability exactly proportional to the degree or strength 
of the target node. In reality, however, this absolute 
quantity information of an agent is often unknown, while 
it is quite common to have a clear idea about the relative 
values of two agents. For this perspective, Fortunato et 
al. propose a criterion of network growth that explicitly 
relies on the ranking of the nodes according to the pres- 
tige measure [9]. This rank-based model can well mimic 
the reality in many real cases that the relative values of 
agents is easier to access than their absolute values. Mo- 
tivated by their work, we propose a model for weighted 
network evolution with only ranking information avail- 
able. Analytically and by simulations, we demonstrate 
that the generated networks possess scale-free distribu- 
tions of degree, strength, and weight. The clustering and 
correlation properties of this class of networks are also 
investigated. 

A weighted network is often denoted by a weighted ad- 
jacency matrix with element Wij representing the weight 
on the link connecting node i and j. In the case of undi- 
rected graphs, weights are symmetric Wij — Wji, as we 
will focus on. A natural generalization of connectivity in 
the case of weighted networks is the node strength de- 
fined as Si = X)jGV(i) ' where the sum runs over the 
set V(i) (neighbors of node i). This quantity is a nat- 
ural measure of the importance or centrality of a node 
in the network. As confirmed by measurement, weighted 
complex network not only exhibits a scale-free degree dis- 




FIG. 1: (Color online) Node strength versus node age with 
a = 0.5, a — 1.0, and a = 2.0 from top to the bottom. Inset: 
Log-Log plot of the relation between node strength and node 
rank by age. All the data are averaged over 100 independent 
runs of network size N = 5000. 



tribution P{k) - fc"'' with 2 ^ 7 ^ 3 [10,11], but also 



the pow-law strength distribution P{s) 



[11] and 



weight distribution P{w) ^ w^" [12]. Highly correlated 
with the degree, the strength usually displays scale-free 
property s ^ [13,14]. 

In the present model, the prestige ranking criterion 
is strength. The definition of the model is based on 
two coupled mechanisms: the topological growth and the 
weights' dynamics. The model dynamics starts from an 
initial seed of No nodes connected by links with assigned 
weight Wo- 

(1) Topological growth. At each time step, a new node 
n is created and m new links, with an assigned weight 
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FIG. 2: (Color online) The inset shows the cumulative 
strength distributions of networks generated by using our 
model with different parameters a = 0.5, a = 1.0, a = 1.5, 
and a = 2.0. The four dashed lines have slopes 2.00, 1.00, 
0.67, and 0.50 separately for comparisons. The main plot 
shows the value of the strength distribution exponent as a 
function of a obtained from numerical simulations. The dot- 
ted line is the prediction of Eq. (6). All the data are averaged 
over 100 independent runs of network size A'^ = 10^. 



these considerations, we check this correlation by numer- 
ical simulations. Fig. 1 shows the node strength as a 
function of its age. It can be found that the function 
is monotone increasing with certain fluctuations. There- 
fore, in the following theoretical analyse, we make an 
approximation that we use the ranking by age instead of 
that by strength. It will be showed by numerical simula- 
tions that this approximation is reasonable. 

The network growth starts from an initial seed of A^o 
nodes, and continues with the addition of one node per 
unit time, until a size A'^ is reached. Hence, each node 
is labeled with respect to the time step of its generation, 
and the natural time scale of the model dynamics is the 
network size N . If the nodes are sorted by age, from the 
oldest to the newest, the label of each node coincides with 
its rank, i. e., Ri = Hi. Therefore, the node strength Si 
is updated according to this evolution equation: 



dsj 
dt 
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Using the continuous approximation, we treat s, w, and 
time t as continuous variables and approximate the sums 
with integrals. Solving Eq. (3) yields the strength evo- 
lution equation: 



Wo to each, are set between node n and pre-existing 
nodes. The previous nodes are ranked according to their 
strength, and the linking probability that the new node 
be connected to node i depends on the rank Ri of i: 



n?! — n 



R- 



(1) 



where a > is a real-valued parameter. Note that the 
larger the rank of the node is, the more difficult for it to 
gain new links, which is reasonable in real life. 

(2) weights' dynamics. Analogous to the step in the 
model proposed by Barrat et al. (BBV model) [15], the 
introduction of the new link on node i will trigger local 
rearrangements of weights on the existing neighbors j £ 
V(i), according to the rule 



Wi 



Wij + d — - , 

Si 



(2) 



where S is the total induced increase in strength of node 
i. 

We firstly investigate the probability distribution of 
the generated network. Since the strength-based ranking 
of a node can change over time, it is hard to analyze the 
model directly by the ranking of node strength. How- 
ever, for a growing weighted network, there is a strong 
correlation between the age of node and its strength, as 
the older nodes have more chances to receive links. For 



(4) 



Consequently, we can easily obtain in the infinite size 
limit the probability distribution: 



Pis) 



(5) 



which shows that the strength distribution of the network 
follows a power law with exponent rj = 1 + 1/a for any 
value of a. 

Similarly to the previous quantities, it is possible to 
obtain analytical expressions for the evolution of weights 
and the relative statistical distribution. The rate equa- 
tion of weight Wij can be written as: 



dwij 
dt 



I — - 



R~ 



Si 



-0 — - + m 



(6) 



Incorporating with Eq. (4), the above equation can be 
solved that Wij {t/tijY^^^~°'\ where = max{i,j) is 
the time at which the edge is established. Therefore, the 
probability distribution P{w) is in this case also a power 
law P{w) ^ , where 
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FIG. 3: (Color online) Insets: Cumulative weight distribu- 
tions of networks built according to the present model for 
different value of (a) S and (b) a. The main plots show the 
value of the weight distribution exponent 6^ as a function of 5 
and a obtained from numerical simulations. The dotted lines 
are the prediction of Eq. (8). All the data are averaged over 
100 independent runs of network size = 10*. 



In order to check the analytical predictions, we per- 
formed numerical simulations of networks generated by 
the present model, where the prestige ranking criterion 
is strength. In the inset of Fig. 2, we plot the cumulative 
strength distributions of the networks corresponding to 
various values of the exponent a. In the logarithmic scale 
of the plot, they exhibit power-law behaviors in agree- 
ment with theoretical results. The relation between a 
and the exponent r/ of the strength distribution is showed 
in the main plot of Fig. 2, which confirms the validity of 
Eq. (6). Together, the power-law distribution of weight 
P{'w) is shown in Fig. 3. The analytical predictions can 
be perfectly confirmed by numerical simulations. Not- 
ing the weights' dynamics step in the definition of the 
model, the triggered increase S is only arranged locally. 
Therefore, we expect the proportionality relation s ^ k, 
by which we easily obtain the scale-free distribution of 
degree P{k) ~ k~'^ with 'y = ij = 1 + 1/a. Since there 
exist no new properties, we do not show them again here. 

To better understand the topology of our model net- 
works, we also studied the clustering and correlation de- 
pending on the model parameters a. The clustering of a 
node i is defined as [3] 



2Ei 



(8) 



where Ei is the number of links between neighbors of 
node i. It measures the local cohesiveness of the net- 
work in the neighborhood of the node. The average over 
all nodes gives the network clustering coefRcient C, which 
describes the statistics of the density of connected triples. 
Further information can be gathered by inspecting the 
average clustering coefficient C(fc), which denotes the ex- 
pected clustering coefficient of a node with k degrees. In 
many networks, the average clustering coefficient C{k) 
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FIG. 4: (Color online) Illustration of the average clustering 
coefficient C as a function of networks size N for (a) a = 
0.5 and (b) a = 2.0. The insets show the behavior of C(fe) 
depending on degree k. The curves in (a) is fit to algebraic 
decay form, 2.50 x A''"" ''^. The solid line in the inset of (b) 
has slope -1 for comparisons. All the data are averaged over 
100 independent runs. 



exhibits a highly nontrivial behavior with a power-law 
decay as a function of k characterizing the intrinsic hi- 
erarchy of the topology [16]. For a = 0.5, the clustering 
coefficient C seems to converge to zero. This is seen by 
the accurate fits to algebraic decay forms in Fig. 4 (a). 
Meanwhile, C(fc) is uncorrelated with fc, denoting that 
the network does not possess hierarchical structure. For 
a — 2.0, C approaches a stationary value of about 0.9 
in the limit of large N , which is showed in Fig. 4 (b). 
In this simple scaling form of clustering coeffi- 

cient, C(fc) ~ is obtained, which indicates that the 
network topology exhibits hierarchical manner. 

Another commonly studied network property is the de- 
gree correlation of node i and its neighbor. The aver- 
age nearest neighbor degree of node with connectivity fc, 
knn{k), is proposed to measure these correlations. If de- 
grees of the neighboring nodes are uncorrelated, knn{k) 
is a constant. When correlations are present, two main 
classes of possible correlations have been identified: as- 
sortative behavior if fc„„(fc) increases with fc, which in- 
dicates that large degree nodes are preferentially con- 
nected with other large degree nodes, and disassortative 
if knn{k) decreases with fc, which denotes that links are 
more easily built between large degree nodes and small 
ones. A simpler measure to quantify this structure is 
assortative mixing coefficient [17]: 



(9) 



where ji,ki are the degrees of nodes at the ends of the 
ith edges, with i = 1,...,L {L is the total number of 
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FIG. 5: (Color online) Illustration of the assortative mixing 
coefficient r as a function of networks size A'' for (a) a = 0.5 
and (b) a — 2.0. The insets show the behavior of fc„„(fc) 
depending on degree k. The curves in (a) is fit to algebraic 
decay form, -0.40 x A^""'^^. The solid line in the inset of (b) 
has slope -1 for comparisons. All the data are averaged over 
100 independent runs. 

edges in the graph). This quantity takes values in the 
interval [—1,1], where positive values mean assortative 
and negative values mean disassortative. Fig. 5 shows 
the simulation results. When a = 0.5, the value of r 
converges algebraically to zero, and knn{k) is unrelated 
with fc, which denotes that correlations are absent. On 



the contrary, when a = 2.0 the assortative mixing coeffi- 
cient is almost independent of network size for large N. 
Meanwhile, fc„„(fc) ^ k~^, characterizing the disassorta- 
tive degree correlation in the network. 

To sum up, we studied a rank-based model for weighted 
network. The scale-free properties of probability distri- 
butions of degree, strength, and weight are obtained ana- 
lytically and by simulation. Furthermore, we investigate 
the clustering and correlation of the network. It is in- 
dicated that a structural phase transition occurs when 
the growth dynamics parameter a = 1 [18]. For a < 1 
(7 > 2), C{k) and knn{k) are observed as a horizontal 
line subject to fluctuations, and clustering coefficient C 
and assortative mixing coefficient r converge to zero in 
the large limit of network size N. For a>l(l<7<2), 
there emerge a few hub nodes in the network which are 
linked to almost every other site ,and the generated net- 
work exhibits hierarchical topology and disassortative de- 
gree correlation. Moreover, the clustering coefficient C 
is independent of network size and approaches a high 
value. Interestingly and specially, in this region of pa- 
rameter a, the generated networks can well mimic the 
biological networks which always appear to be disassor- 
tative [3,17] and possess hierarchical organization [16,19]. 
We think that this class of network provides us with a 
new method to reconstruct the hierarchies and organiza- 
tional architecture of biological networks, and it may be 
beneficial for future understanding or characterizing the 
biological networks. 
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